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ABSTRACT 



6 

+3 • We have studied three-dimensional hydrodynamic interactions of relativistic 

extragalactic jets with two-phase ambient media. These jets propagate through 
a denser homogeneous gas and then impact clouds with densities 100 to 1000 
times higher than the initial beam density. The deflection angle of the jet is 
influenced more by the density contrast of the cloud than by the beam Mach 
number of the jet. A relativistic jet with low relativistic beam Mach number can 
eventually be slightly bent after it crosses the dense cloud; however, we have not 
seen permanently bent structures in the interaction of a high relativistic beam 
Mach number jet with a cloud. The relativistic jet impacts on dense clouds do 
not necessarily destroy the clouds completely, and much of the cloud body can 
survive as a coherent blob. This enhancement of cloud durability is partly due to 
the geometric influence of the off-axis collisions we consider and also arises from 
the lower rate of cloud fragmentation through the Kelvin-Helmholtz instability 
for relativistic jets. To compare our simulations with observed extragalactic ra- 
dio jets, we have computed the approximate surface distributions of synchrotron 
emission at different viewing angles. These surface intensity maps show relativis- 
tic jets interacting with clouds can produce synchrotron emission knots similar 
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to structures observed in many VLBI-scale radio sources. We find that the syn- 
chrotron emission increases steeply at the moment of impact and the emission 
peaks right before the jet passes through the cloud. 

Subject headings: galaxies: active — galaxies: jets — hydrodynamics — ISM: 
clouds — methods: numerical — relativity 

1. Introduction 

Relativistic jets emerging from extragalactic sources associated with active galactic nu- 
clei (AGNs) are the most important means of transporting energy and mass from AGNs to 
an external medium over large distances. To understand how these relativistic jets interact 
with an inhomogeneous external medium containing small, dense gas clouds or clumps has 
been recognized as important for a long time. These interactions may change substantially 
the direction of relativistic jet flows, trigger extensive star formation in the shocked clouds, 
and possibly explain the basic mechanism behind the morphology of many extragalactic 
radio jets (Fanaroff & Riley 1974). 

Recent observations have revealed strong evidence of features associated with changes 
in jet directions resulting from interactions with small gas clouds in the narrow-line regions 
of Seyfert galaxies (e.g., Mundell et al. 2003). Fast outflows of gas observed in the central 
regions of powerful radio galaxies can also be caused by such interactions (e.g., Emonts et 
al. 2005; Morganti et al. 2005). The most likely interpretation of fast outflows is that all 
gas clouds are not destroyed by the jet; some clouds can severely disrupt the jet while some 
clouds are accelerated to the observed high outflow velocities by the thrust of the jet. It is 
argued that despite of the high energies involved in the interactions, only a few percent of 
the outflowing gas appears to be ionized, while the rest of the gas cools and becomes neutral 
due to highly efficient cooling near the jet bow shock. 

In the context of nonrelativistic hydrodynamic simulations, previous numerical works 
were performed to investigate jet interactions with clouds (de Gouveia Dal Pino 1999; Hig- 
gins et al. 1999; Wang et al. 2000; Saxton et al. 2005) or jets crossing a medium interface 
(e.g., Wiita et al. 1990; Wiita & Norman 1992), focusing on the effects of the interactions on 
the morphology and kinematics of jets. Others studied shock interactions, focusing on the 
structure and evolution of the clouds produced by the interactions in adiabatic cases (Klein 
et al. 1994; Xu & Stone 1995; Poludnenko et al. 2002) and in radiative cases (Mellema et 
al. 2002; Fragile et al. 2004). According to de Gouveia Dal Pino (1999), simulations with 
conditions appropriate to protostellar jets making off-axis collisions with clouds produced 
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a deflected beam. The deflection angle tended to decrease with time as the beam slowly 
penetrated the cloud and when the jet penetrated most of the cloud, the deflected beam 
faded and the jet resumed its original propagation direction. Wang et al. (2000) found the 
following: powerful extragalactic jets eventually destroyed the clouds they considered, and 
these collisions induced nonaxisymmetric instabilities in the jets; weak jets can be effec- 
tively halted or destroyed by massive clouds; and slow, dense jets that were bent remained 
stable for extended times. Synthetic radio images produced by hydrodynamic simulations 
for comparison with observations also supported the hypothesis that these interactions are 
responsible for the distorted structures of some radio jets (e.g., Higgins et al. 1999). All 
those numerical works considered nonrelativistic jet speeds less than 0.5c, but the observed 
apparent superluminal motions of extragalactic radio sources indicate intrinsic jet speeds up 
to at least 0.98c (Zensus 1997). Thus, it is essential to perform relativistic hydrodynamic 
simulations of this problem in order to cover the range of true jet speeds. 

Since time-dependent numerical simulations of relativistic jets were first reported (van 
Putten 1993; Duncan & Hughes 1994; Marti et al. 1994), multidimensional relativistic hydro- 
dynamic simulations have been used as an important method in understanding of relativistic 
jets (Marti et al. 1997; Komissarov & Falle 1998; Aloy et al. 1999; Rosen et al. 1999; Hughes 
et al. 2002; Mizuta et al. 2004). The morphological and dynamical properties of relativistic 
jets propagating through a homogeneous medium were studied by Marti et al. (1997) in two 
dimensions and by Aloy et al. (1999) in three dimensions. Komissarov & Falle (1998) inves- 
tigated the large-scale flows produced by classical and relativistic jets in a uniform external 
medium using analytical and numerical studies. Hughes et al. (2002) performed in three 
dimensions a study of the deflection of relativistic jets by an oblique density gradient and of 
the precession of relativistic jets. They found that fast relativistic jets can be significantly 
influenced by an oblique density gradient, showing a rotation of the Mach disk with the flow 
bent via a strong oblique internal shock. 

In this paper we present results from three-dimensional hydrodynamic simulations of 
the interactions of relativistic jets with dense clouds. We focus on the off-axis collision of 
the relativistic jet with a steady spherical cloud. The main concerns of this study are how 
the relativistic jets are influenced by these interactions and how the interaction affects the 
evolution of the cloud. In §2 we outline briefly the dynamical problem, while the basic 
equations, numerical method and setup we employ are described in §3. In §4 we describe 
the results, and we present a summary and discussion in §5. 
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2. Problem Description 

We consider the three-dimensional interactions of relativistic jets with two-phase ambi- 
ent media. These jets propagate through a denser ambient gas and then hit spherical clouds 
with densities higher than that of the ambient gas. The initial ratio of the cloud density, p c , 
to the ambient medium density, p a , and that of the beam density, p b , to the ambient medium 
density, are respectively defined as 

_ Pc _ Pb m 

Pa Pa 

If we neglect complicating effects including radiative cooling and gravity and consider only 
hydrodynamic effects, then this problem can be relatively simple and depends only on a few 
hydrodynamic variables, the Mach number of the jet and the initial density contrasts given 
in equation (1). Any geometric effects, such as different impact zone sizes or cloud shapes, 
certainly will make differences in the evolutions of jets and clouds, but the overall dynamical 
evolutions should not be very sensitive to them. Thus, we will focus on the evolutions of jets 
and clouds influenced by the above hydrodynamic effects. 

The approximate propagation velocity of the jet through the homogeneous ambient 
medium can be obtained by the conservation of the momentum flux of the beam and ambient 
medium in the reference frame of the Mach disk (Marti et al. 1997). Assuming pressure 
equilibrium between the beam and the ambient medium, the conservation of the momentum 
flux is p b h b T' b 2 v' b 2 = p a h a T' 2 v' 2 with the following relations, v' b = (v b - v a )/(l - v b v a ), T' b = 
T b T a (l — v b v a ), v' a = —v a , and T' a = T a . Here h is the specific enthalpy and v' and V 
represent respectively velocity and Lorentz factor measured in the reference frame of the 
Mach disk, while v and T indicate those measured in the rest frame of the ambient medium. 
The subscripts b and a stand for the beam and the ambient medium, respectively. After 
substituting for the primed variables in terms of the unprimed ones, the conservation of the 
momentum flux is derived to be 

p b h b T 2 b (v b - v a f = p a h a v 2 a . (2) 

Then the one-dimensional jet advance velocity, estimated in the rest frame of the ambient 
medium is 



i?*5r (4) 

Pa^a 



where r/* is given by 

V 

In the nonrelativistic limit (h — > 1, T — > 1), rj* approaches i], so that v a represents the 
classical jet advance velocity through the ambient medium, i.e., v a = v b /(^l/r] + 1). 
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Based on this jet advance velocity, we define the dynamical timescale called the beam 
crossing time t hc , 

the = — , (5) 

v a 

as the time taken for the beam to sweep a distance across the ambient medium equal to 
the diameter of a cloud with radius r c . Since the timescale t bc depends only on a single 
variable v a (for fixed cloud radius), it is extremely useful in comparing and characterizing 
the dynamical evolutions of both jets and clouds with different model parameters. 

Although we use the beam crossing time as the primary timescale in this study, it is 
also interesting to estimate the cloud crushing and cooling timescales. The cloud crushing 
timescale is the time required for the beam to cross the cloud diameter during the phase 
of cloud compression, and if v a is nonrelativistic, this timescale can be approximated as 
t cc ~ 2^/x~r c /v a (Klein et al. 1994). Clearly, t bc ~ t cc in the absence of clouds, and for dense 
clouds (x ^> 1), tbc < t cc . Following Fragile et al. (2004), the cloud cooling timescale can be 
roughly estimated from t coo i ~ Cv^/(x 3 ^ 2 p c ), where the constant C = 7.0 x 10~ 35 g cm -6 s 4 . 
With values reasonable for kpc-scale extragalactic situations, r c = 1 kpc, v a = 0.1c, x = 100, 
and p c = 10 2 m H cm -3 , we find that tbc < t cc ~ t coo i. Thus cooling will not be extremely 
important during the cloud compression phase for the chosen values. For fixed density and 
cloud radius, the cloud cooling timescale becomes longer compared to the cloud crushing 
timescale as v a increases, so the effect of cooling is somewhat reduced for relativistic jets 
compared to nonrelativistic ones. For parameters more relevant to VLBI-scale jet/cloud 
collisions, r c = 0.5 pc, v a = 0.5c, x — 10 4 , and p c = 10 6 m H cm~ 3 , we have tbc ~ t coo i < t cc , 
so cooling would be more important in this case. A more detailed consideration of cooling 
timescales is beyond the scope of the current paper. 

Three distinct evolutionary stages can be considered in this problem. There is an initial 
jet propagation stage where the jet advances through a homogeneous ambient medium with 
velocity v a . Once a jet is launched, a bow shock propagates into the ambient medium; 
this is followed by a Mach disk shock in the beam which is quickly established during this 
stage. When the jet strikes the cloud, the jet transmits a shock into the cloud. If the 
cloud/ambient density contrast is sufficiently large and the jet speed is relatively slow, the 
speed of the transmitted shock in the cloud is much slower than that of the bow shock of the 
jet. Thus the bow shock entirely encloses the cloud, which leads to the development of the 
Kelvin- Helmholtz instability at the cloud surface (e.g., Klein et al. 1994). The final stage 
is when the jet passes through the cloud. In this phase the cloud begins to reexpand just 
after the jet reaches the rear edge of the cloud. At the same time, the jet propagates in the 
original direction or in a new direction. 
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3. Numerical Simulations 

3.1. Basic Equations 

The special relativistic hydrodynamic equations are written in a covariant form (e.g., 
Landau & Lifshitz 1959; Wilson & Mathews 2003) 

d a ( P U a ) = 0, (6) 

d a T a<i = 0, (7) 
where the energy momentum tensor is given by 

T «p = ( e + p ) u a l/P + pg a/3 . (8) 

Here, d a = d/dx a is the covariant derivative with spacetime coordinates x a = [t,Xj], U a = 
[T,Tvj] is the normalized (U a U a = — 1) four-velocity vector, and a metric tensor g a ^ with 
a signature +2 is used. The mass density, velocity, internal plus mass energy density, and 
pressure in the local rest frame are denoted by p, Vj, e, and p, respectively. Greek indices 
(e.g., a, (3) denote the spacetime components while Latin indices (e.g., i, j) indicate the 
spatial components, and the speed of light is set to unity (c = 1) throughout this paper. 

For our numerical purposes, it is convenient to rewrite the covariant equations (6)-(8) 
in the index form which gives a hyperbolic system of conservation equations 

dMi d , 

-df + dx~ {MtV > +p6l ' ) = °> (10) 

_ + _[(* + „)„,] = 0, (ID 

where the equation of state (EOS) is given by 

p = ( 7 - 1) (e - p) . (12) 

Here, D, Mj, and E are the mass density, momentum density, and total energy density in 
the reference frame, respectively, and 7 is the adiabatic index. We note that we restrict 
ourselves to an ideal gas EOS in this study (cf. Ryu et al. 2006). 

The quantities in the reference frame are related to those in the local rest frame via 
following transformations 

D = Tp, (13) 
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M t = T 2 (e + p)v l , (14) 
E = T 2 (e + p) - p, (15) 

where the Lorentz factor is given by 

1 



(16) 



with v 2 = v\ + v 2 + v\. 



If an EOS is assumed, the local sound speed c s and the specific enthalpy h are easily 
derived. For an ideal gas, they are given by 



1 dp i dp u n IP 



c = i? + ^> fc = i + — V 1 - (it) 

hop oe 7 — 1 p 



A 7- law gas such as an ideal gas has the local sound speed limit c s < y^y — 1 and in the 
ultrarelativistic case e ^> p, the local sound speed approaches its limit (i.e., c s — > ^7 — 1). 



3.2. Numerical Method and Setup 

The system of equations (9)-(12) can be solved numerically with explicit finite difference 
upwind schemes which are based on exact or approximate Riemann solvers using the char- 
acteristic decomposition of relativistic hydrodynamic conservation equations. Although the 
upwind schemes were originally developed for nonrelativistic hydrodynamics, some schemes 
have been extended to special relativistic hydrodynamics retaining the advantages of the 
upwind schemes, including high accuracy and robustness. 

A multidimensional code for solving the special relativistic hydrodynamic equations as 
a hyperbolic system of conservation laws based on the total variation diminishing (TVD) 
scheme (Harten 1983) was developed and tested in Choi & Ryu (2005). The TVD scheme is 
an explicit Eulerian finite difference upwind scheme and an extension of the Roe scheme to 
second-order accuracy in space and time. The code uses a new set of conserved quantities, 
which lead to a new eigenstructure for special relativistic hydrodynamics and employs an 
analytic formula for the calculation of the local rest frame quantities from the reference frame 
quantities. The advantage of our code is that it is simple and fast, and yet it is accurate and 
reliable enough. The performance of the code was demonstrated through several standard 
tests, including relativistic shock tubes, a relativistic wall shock, and a relativistic blast 
wave, as well as test simulations of the relativistic version of the Hawley-Zabusky shock 
and a relativistic extragalactic jet (Choi & Ryu 2005). For our new simulations, we have 
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parallelized this code using the message-passing interface (MPI). The simulations described 
here typically use 64 processors on a Linux cluster. 

Table 1 lists the initial parameters of the four different relativistic jet-cloud interaction 
models we have investigated in this study. All models use the adiabatic index 7 = 4/3 and 
assume pressure-matched jets and clouds, i.e., Pb/p a = Pc/p a = 1, where pb, p c , and p a are 
the pressure of the beam, cloud, and ambient medium, respectively. We set c = r c = p a = 1 
in our models, so that all physical quantities are dimensionless and can be scaled to any 
specific physical model (e.g., t — > tc/r c , p — > p/p a )- The Newtonian beam Mach number, 
Ai^ = Vb/c Sy b, where c s ^ is the sound speed in the beam, as well as the relativistic beam 
Mach number, M.f = {T b /^ s,b)-M-b > where T s ^ is the Lorentz factor associated with the 
beam sound speed are given in Table I. As discussed in Konigl (1980), in the context of 
relativistic gasdynamics the relativistic Mach number is the best analog of the Newtonian one 
for nonrelativistic flows, so we use this relativistic beam Mach number to describe physical 
properties in our models. The initial density contrast between the beam and the ambient 
medium is fixed to rj = 0.1, so that jets strike clouds with densities 100 to 1000 times higher 
than the initial beam density. In models Ml and M2, the clouds interact with the lower 
relativistic beam Mach number jets, where the relativistic effects are less dominant, with 
smaller beam velocities and internal energies. Model Ml is almost identical to model M2 
except for the smaller density ratio of the cloud to the ambient medium. Models M3 and 
M4 have been chosen to study the cloud interactions with the higher relativistic beam Mach 
number jets, with the more dominant relativistic effects caused by larger beam velocities and 
internal energies. Again, the initial conditions of model M3 are the same as those of model 
M4 except for the smaller density ratio of the cloud to the ambient medium. 

We set up the density gradient of the spherical cloud edge with a hyperbolic tangent 
function 



where r is the distance from the center of the cloud and Ar is the scale parameter for the 
width of density transition (Ar <C r c ). The presence of a true density discontinuity instead of 
this steep function would not affect the dynamics of jet-cloud interactions significantly, but 
the discontinuous cloud edge is approximated by this somewhat smoothed density profile to 
avoid numerical artifacts as the jet impacts the cloud. We assume other physical quantities 
such as pressure and velocity are constant across the transition width. 

The simulations have been performed in the three-dimensional computational domain 
with x = [0,8], y = [0,8], and z = [0,8] using a uniform Cartesian grid of 256 3 cells. 
The beam, with a circular cross section of radius r b = 1/4 (8 cells), is initially located 
at (x, y, z) = (0, 4, 4) and propagates through the ambient medium along the positive x- 




(18) 
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direction. In order for the relativistic jet to collide off axis with the cloud at rest, the center 
of the cloud, with radius r c — 1 (32 cells), is placed at (x,y,z) = (4,3.5,4); hence the 
relativistic jet hits the spherical cloud with an impact angle of 30°. The outflow boundary 
condition is imposed on all boundaries of the computational domain except where the inflow 
boundary condition is used to maintain the continuous jet. We were able to assign the 
relativistic jet only 8 cells per initial beam radius and the cloud 32 cells per initial cloud 
radius due to the limitation of computational resources. This resolution is less than that 
of previously reported two-dimensional works which are related to this problem. Thus our 
three-dimensional simulations may not be fully converged and some quantities to be described 
may change if three-dimensional simulations with much higher resolutions are performed in 
future studies; however, our tests of the code do indicate that simulations with this level of 
resolution should be reasonably accurate (Choi & Ryu 2005). 



4. Results 

4.1. Morphology and Dynamics 

The gray-scale images in Figures l(a)-(d) show the distinct evolutionary phases of 
models M1-M4, respectively. These images show the x — y plane with z = 4 in the three- 
dimensional computational domain. In each of these figures the top to bottom panels rep- 
resent density, pressure, and Lorentz factor, respectively (in logarithmic scales) while the 
left to right panels represent evolutionary stages shown at three different times, t = t^ c , 
(the + ^cnd)/2, and t cnd . 

The early stages of the relativistic jet propagation through the uniform ambient medium 
until the jet is about to collide the cloud (t/t hc ~ I) are basically similar to those found 
in earlier simulations (e.g., Marti et al. f997; Aloy et al. 1999). Several key features are 
clear from the left panels of Figures l(a)-(d). In all the models a bow shock that separates 
the jet from the external medium is driven, the beam itself is terminated by a Mach disk 
(terminal shock) where the beam kinetic energy is converted into its internal energy, and 
shocked jet material flows backward along the contact discontinuity (working surface) into a 
cocoon. There is no difference between models Ml and M2 and between models M3 and M4 
at this stage because of the same initial conditions of the jets and the same ambient media 
properties for these two pairs of simulations. 

The relativistic beam Mach number of the jet is associated with the shape of the bow 
shock. In models Ml and M2, the lower relativistic beam Mach number jets, with a lower 
propagation velocity (v a ~ 0.42) and internal energy, have bow shocks with narrower conical 
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shapes, and the Mach disk is quite close to the bow shock. This conical shape of the bow 
shock tends to be broader as the relativistic beam Mach number of the jet increases, as 
seen for models M3 and M4; these higher relativistic beam Mach number jets, with a higher 
propagation velocity (v a ~ 0.78) and internal energy, also have the Mach disk standing off 
farther from the bow shock. The shapes of the bow shocks are also connected with the sizes 
of the impact cross section when the jets begin to interact with the cloud. The low relativistic 
beam Mach number jets in models Ml and M2 feature relatively thick cocoons while the 
high relativistic beam Mach number jets in models M3 and M4 have thin cocoons. This 
dependence of the cocoon morphology on the relativistic beam Mach number is consistent 
with previous results (see e.g., Marti et al. 1997). Although the structural differences in 
the jet head and the cocoon are evident by this early stage of the evolution, the internal 
structures within the beam and backflows are not dominant and are barely distinguishable 
at this stage. 

In every model, the relativistic jet begins to partially deflect as a direct response to its 
interactions with the clouds. This is seen in the middle column of panels of Figures l(a)-(d); 
seen most clearly in the middle bottom panels are fast streams coming from the Mach disk 
at significant angles with respect to the jet axis. These deflection features are stronger in 
models M2 and M4, which have higher ratio of the density of the cloud to that of the ambient 
medium (x — 100). The deflection angles of the portion of the post-Mach shock flows with 
respect to the beam propagation axis are very time-dependent. In our models these angles 
peak when the jets cross over approximately half the clouds (at t/t hc ~ 2.5, 3.5, 2.5, and 
3 for models M1-M4, respectively). For the comparable dynamical times, models M2 and 
M4, both with \ = 100 but having different beam Mach numbers, show 80° — 90° deflection 
angles, while models Ml and M3 with the same beam Mach numbers as the models M2 
and M4, respectively, but with x — 10, show smaller deflection angles of about 45°. This 
indicates that the deflection angle is more strongly influenced by the density contrast, Xi 
than by the beam Mach number of the jet. For an off-axis collision there are weak deflection 
features on the other side of the jet axis, where the deflection of the outflow from the beam 
is significantly suppressed by the dense cloud. That suppression leads to the production 
of a strong oblique shock within the beam. As seen in the figures, the oblique shocks are 
quite strong in models M2 and M4, but in models Ml and M3 there are only relatively weak 
oblique shocks in the beam. Comparing at this stage models Ml and M3 with models M2 
and M4, we note that the bow shocks enclose less of the cloud in models Ml and M3 because 
of their lower density contrast, x- That implies quicker penetration of the clouds by these 
jets, so the strengths of the oblique shocks in these beams are reduced. 

Some additional properties of the simulations at this stage are shown in Figures 2 and 3. 
Figure 2 illustrates one-dimensional flow structures of density, pressure, and Lorentz factor 
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along the beam propagation axis for models Ml and M3 at the same epoch as in Figures 
1(a) and (c). In both models there are spikes in the density and pressure associated with 
the impact by the incident jets, while there is little change in beam Lorentz factor. Figure 3 
shows the images of the logarithm of the Lorentz factor projected at the viewing angle of 0° 
for models M2 and M4, at t/tbc = 3.5 and 3, respectively. These projection images clearly 
show the anisotropic distribution and directions of the deflected gas induced by the jet. This 
gas is an admixture of jet and cloud material, but only a small fraction of the cloud gas is 
shown in these projection images since the mean cloud velocity computed in each component 
(refer to §4.2) is (vi) < 0.01 and 0.06 for models M2 and M4, respectively, at the same epoch 
as in Figure 3. This implies that this deflected gas consists predominantly of jet material 
although small amount of cloud material is entrained in these deflected structures. This 
presence of deflected gas accelerating toward a terminal velocity strongly suggests that such 
deflected and accelerated gas is responsible for at least some of the outflowing gas observed 
in the vicinity of AGNs. 

Once the jet passes through the cloud, it begins to accelerate, causing a change in the 
shape of the bow shock. As visible in the right panels of Figures 1 (a)-(d) shown when the 
jet head nearly reaches the boundary of the computational cube (at t/% c = 4, 6, 4, and 5 for 
models M1-M4, respectively), the shape of the bow shocks changes more clearly in the low 
relativistic beam Mach number jets than in the high relativistic beam Mach number jets. 
That reflects the fact that the acceleration of the jets is somewhat faster in low relativistic 
beam Mach number jets. That reacceleration occurs in essentially the original propagation 
direction or in a somewhat new direction. In our simulations there is a trend for the flow 
of the jet to be bent more when a lower relativistic beam Mach number jet interacts with a 
denser cloud, with the least bending seen for model M3 and the most for model M2. We see 
in the right panels of Figure 1(b) that the beam is bent by about 10° with respect to the 
original jet axis. The bent jet still remains stable and collimated over the several dynamical 
times we could follow its development. 

After the jet head passes the cloud, the amount of strongly deflected gas gradually 
reduces and the oblique shocks continue to develop in the beam. These oblique shocks are 
unlikely to play a major role in slowing the jets because we do not find any significant 
deceleration features during this stage. Although a significant portion of the momentum 
flux of the jets is transfered to the deflected gas and the cloud through the collision events, 
the jets in our simulations are still stable and well collimated over several dynamical times 
after collisions even if the jet is bent. This stable, collimated condition is quantitatively 
apparent in the flow structures of density, pressure, and Lorentz factor shown in Figure 2 
at t/tbc = 4 (corresponding to the dashed lines) for models Ml and M3, respectively. There 
are only slight fluctuations in the flow structures at this late stage. 
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In comparing our simulations with hydrodynamic simulations of nonrelativistic jet- 
cloud interactions (de Gouveia Dal Pino 1999; Higgins et al. 1999; Wang et al. 2000) we 
can only note some fairly basic similarities and differences between our relativistic models 
and the roughly corresponding nonrelativistic models. The lack of good overlap between the 
M.^ values for the relativistic jets and the standard Mach number for the nonrelativistic 
jets as well as differences between cloud size to jet-width ratios considered here and in that 
earlier nonrelativistic work prevents us from making quantitative comparisons. Relativistic 
jets interacting with dense clouds certainly do show general morphological features such as 
deflections of some gas and bent structures of jets similar to those found in some of the 
nonrelativistic jet-cloud interactions. The slower relativistic jet shows a bent structure after 
interaction, which is similar to that found in nonrelativistic simulations involving "weak" 
jets while the faster relativistic jet effectively plows through the clouds. Higher power non- 
relativistic jets also can plow through, and apparently completely destroy, clouds. However, 
some major differences arise because of the larger propagation velocity of the relativistic jets. 
Because of this, moderately light (i] = 0.1) relativistic jets are not effectively decelerated and 
disrupted by the dense (x = 100) clouds, whereas nonrelativistic jets assaulting clouds of 
similar density ratios typically are disrupted. Our relativistic jets are rather reaccelerated in 
either a slightly new or essentially the original direction after their interactions with clouds. 
The large propagation velocity also suppresses the development of hydrodynamic instabili- 
ties in the jets, so that the jets still remain stable and collimated even after the jets smash 
into much denser clouds. In addition, as discussed in §4.2, clouds impaled by relativistic jets 
also appear to survive somewhat better than do those hit by strong nonrelativistic jets. 



4.2. Cloud Evolutions 

Although previous studies of jet interactions with clouds mainly emphasized the dynam- 
ical and morphological features of the jet itself, it is also important to follow the evolution 
of the clouds during and after the off-axis collisions with relativistic jets. One key reason for 
investigating the fate of the clouds is that the leftover cloud material is a strong candidate 
for star formation regions in the vicinity of AGNs (e.g., Rees 1989; Gopal-Krishna & Wiita 
2001). The cloud is expected to undergo a somewhat different evolution in our case com- 
pared with the evolution of the clouds struck by the nonrelativistic planar shocks considered 
in earlier work (e.g., Klein et al. 1994; Xu & Stone 1995; Mellema et al. 2002; Fragile et al. 
2004). 

In order to describe the evolution of a cloud quantitatively, we introduce a conserved 
variable / called a Lagrangian tracer (e.g., Jones et al. 1996) which is updated along with 



the primitive hydrodynamic variables in our simulations. The evolution of the Lagrangian 
tracer is followed by the conservation equation 



dDf 



+ 







(Df Vj ) = 0. 



(19) 



dt 



Since the above conservation equation is almost identical to the mass conservation equation 
(9), it is separately solved using the same TVD routine adopted for solving the mass conser- 
vation equation. Initially the tracer variable is set to unity (f c = 1) inside the cloud while 
the variable is set to zero (/ c = 0) everywhere outside the cloud, so that the density of cloud 
material is given as D c = Df c (i.e., p c = pf c ) for a given tracer f c in any zone. Then the 
total mass of the cloud is computed by the integration over the entire volume V, 



Jv 

where dV = dxdydz. This enables us to compute the several useful mass-weighted quantities 
such as the mean square radius of the cloud and the mean velocity of the cloud, 



The index % given above represents each spatial component. Another useful mass- weighted 
quantity is the mean thermal energy inside the cloud (e t h)- This is also computed using the 
same volume integration given above. 

We show in Figure 4 the volume-rendering images of cloud density for model M4 at three 
different times, t/tbc = 1, 3, and 5. As a direct consequence of the impact on the cloud by 
the jet the cloud develops a cavity in the cloud body as shown in the figure. The cloud cavity 
continues to grow until the jet completely penetrates the cloud, elongating the cloud material 
outside the cavity along the bow shock of the jet. Unlike the cases studied earlier where a 
cloud interacts with a plane-parallel shock (Klein et al. 1994; Xu & Stone 1995), the cloud 
material is not completely destroyed by the impact of the jet. Some cloud mass is carried 
into the deflected material of the jet, eroding the cloud body, but much of the cloud mass 
remains in a large, coherent blob for at least a few beam crossing times. This enhancement 
of the cloud durability is apparently primarily due to the geometric influence of an off-axis 
collision. Computational resource limits prescribe that we can accurately investigate the 
clouds for only a few beam crossing times, which is less than the many dynamical times for 
which it would be desirable to follow their evolutions. 
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Figure 5 shows for every model the time evolutions of the root mean square radius of 
cloud, the mean cloud velocity, and the mean thermal energy of the cloud. In every model 
the clouds remain in the initial root mean square radius (rj) 1 / 2 = 0.44 until t/% c ~ 1.5. 
When the jet hits the cloud, the cloud is first crushed in the x-direction, along which the jet 
propagates, and then it begins to expand beyond its initial size. The initial compressions in 
the y- and z-directions are very small and the cloud soon gradually expands in both these 
transverse directions. By the end of these simulations the root mean square radii of the 
clouds have expanded to about 1.5 — 2 times their initial values. 

After t/tbc ~ 1-5 the high pressure inside the cloud generated by the incident jet causes 
the entire cloud to accelerate. Unsurprisingly, the acceleration is faster in the x-direction for 
the faster jets in models M3 and M4 and for the lighter clouds in models Ml and M3. The 
mean velocity of the clouds peaks at values between 0.01 and 0.15 in the x-direction and 
0.005 and 0.05 in the y-direction. Note that the mean velocity of the clouds in the z-direction 
is zero because of symmetry in this direction. The maximum velocity of the cloud is always 
rather modest even if the incident jet has a relativistic speed, though if we had considered 
less massive clouds they obviously could have been accelerated to higher speeds. 

As we expect, the mean thermal energy of the cloud increases while the jet strikes the 
cloud. The maximum mean thermal energy of the cloud reaches about 5 — 15 times its initial 
value, depending upon the model. In each model the peaks of the mean thermal energy 
inside the cloud and the mean velocity of the cloud take place nearly at the same time. Note 
that at this point the cloud reexpands after the cloud reaches the maximum compression in 
the ^-direction. 

As mentioned earlier, the jet interaction with the cloud shows that the beam pene- 
trates through the cloud body, which may begin a fragmentation process. A strong shear 
layer developing at the cloud boundary as a result of the interaction with the jet may lead 
to Kelvin-Helmholtz instabilities which enable the disrupted cloud body to fragment. So 
eventually the gas cloud might be broken into small pieces. However, the Kelvin-Helmholtz 
instability becomes inefficient if the density contrast of two slipping fluids is large or if the 
flow is supersonic (Chandrasekhar 1961), so we may not see rapid fragmentation in the 
clouds. Although the fragmentation timescale is difficult to estimate, our simulations show 
no significant cloud fragmentation by t/t^ c ~ 4 — 6. This indicates that the high density 
contrast between clouds and beams and supersonic velocity of the clouds induced by the 
relativistic jets do indeed lower the growth rate of the Kelvin-Helmholtz instabilities. 
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4.3. Synchrotron Emission 

Propagating relativistic jets produce nonthermal radio (synchrotron) emission which 
originates from relativistic high-energy particles accelerated at the shock front. Jones et al. 
(1999) and Tregillis et al. (2001) calculated the synchrotron emission in extragalactic jets by 
explicitly calculating the acceleration of electrons at shocks and following the evolution of 
magnetic field. However, they assumed nonrelativistic jets, and hence the emissivity needs 
to be further examined using relativistic jets. To compute the synchrotron emission from 
relativistic jets, other relativistic hydrodynamic simulations have worked with a simpler 
approximation (Gomez et al. 1997; Komissarov & Falle 1997; Mioduszewski et al. 1997; Aloy 
et al. 2000). Using this same simple model we now calculate the synchrotron emission in 
our simulations in order to estimate how the relativistic jet interaction with a cloud would 
appear in emission as a extragalactic radio source. We make the usual assumptions that the 
jet is optically thin and only the jet material radiates. Thus, in order to separate the jet 
material from the ambient medium and the cloud, we include an additional tracer variable 
fb (see §4.2) which is initially set to unity inside the jet (fb = 1) and zero everywhere outside 
the jet (f b = 0). 

The relativistic high-energy electrons responsible for the synchrotron emission are as- 
sumed to have a power-law energy distribution. Given the spectral index a, the high-energy 
particle number density N , and the magnetic field intensity B, the synchrotron emissivity 
at frequency v is then approximated by the power-law distribution (see e.g., Mioduszewski 
et al. 1997) 

j v oc N B a+1 u~ a . (23) 

The high-energy particle number density, Nq, is assumed to be proportional to the relativistic 
electron energy density, u e , from the integration of the power-law energy distribution over 
some energy range, and u e is also taken to be proportional to the hydrodynamic pressure. 
Then we have N oc u e oc p. Assuming that there is an equipartition of the magnetic field 
energy density ub and the relativistic electron energy density {ub = u e ), then ub oc p. This 
leads to B oc u 1 ^ 2 oc p x l 2 . Therefore, equation (23) becomes 

j v <xp {a+3)/2 v~ a . (24) 

This equation shows that the local thermal pressure approximately reflects the local syn- 
chrotron emissivity. We have used a = 0.75 in our calculation. By integrating the syn- 
chrotron emissivity along the line of sight L at a viewing angle 9, we can compute the 
synchrotron intensity on the surface projected onto the line of sight at the viewing angle 



I u = Jv 2 ]u dL, (25) 
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where the Doppler boosting factor is given by 

V = — -. (26) 

Other relativistic effects, including light aberration and time dilation, have not been included 
in this calculation, as we assume that these effects are negligible. 

Figure 6 shows the synchrotron intensity maps of models Ml and M2 at the viewing 
angles 90°, 45°, and 0°. These maps are shown at t/% c = 2.5 for Ml and 3.5 for M2 when 
the jet is colliding with the cloud. The peak intensity in this figure varies with the models 
and the angles of view. Doppler boosting has a little effect on the emission of the jet at 
the viewing angle 90°, so that the observed emission is very closely related to the intrinsic 
emissivity in this case. At smaller viewing angles (e.g., 45° and 0°), however, the emission 
morphology is determined to a large degree by Doppler boosting. The synchrotron emission 
is dominated by the bright hotspot, which takes the form of the compact emission knot in 
VLBI radio maps. Although the beam and the deflected material show only weak emission 
features, there is a faint secondary spot seen from deflected material in model M2. 

The time evolution of the total synchrotron intensity for models M1-M4 at the viewing 
angles of 90°, 45°, and 0° are shown in Figure 7. The total synchrotron intensity computed 
here is in arbitrary units. There are significant quantitative differences among the models, 
but the intensity curves show qualitatively the same trends. The total synchrotron intensity 
is much amplified at smaller viewing angles of 45° and 0° because Doppler boosting plays 
a role in the amplification of the intensity in these cases. As expected, the passage of the 
jet over a cloud enhances the synchrotron intensity; there are high amplitude bumps in the 
intensity curves during the interactions. The total intensity steeply increases at the moment 
of the impact by the jet, and then gradually increases until the jet crosses over the cloud. 
This tells us that the compression of the plasma in this region produces higher synchrotron 
emission in this approximation where it is tied to the pressure. The peak synchrotron 
intensity occurs shortly after the jet passes through the entire cloud, and after that the 
intensity falls off slowly because the compression is weaker. 

Although we have not computed the thermal X-ray emission in detail, we can briefly 
discuss it. Since the free- free emission (bremsstrahlung) is proportional to p 2 , the total X-ray 
luminosity due to thermal bremsstrahlung is most sensitive to the density of gas, provided 
that the gas is, or becomes, hot enough to emit X-rays. The relativistic jet itself is not 
expected to emit thermal X-rays because of its low density, and only in some cases does the 
synchrotron spectrum extend far enough to produce nonthermal X-ray emission (Harris & 
Krawczynski 2006). The dense cloud is very unlikely to start out hot enough to emit X-rays. 
However, during the jet-cloud interaction the density and pressure of the cloud become so 
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high that the total X-ray emission may be larger inside the cloud than elsewhere, for example, 
in the bow shock of the jet. An estimate of the increase in the total X-ray luminosity of 
the cloud is given by L x /L xfl ~ (pj Pc,o) 2 (T c /T cfi ) 1/2 (y c /V cfi ), where L x is the total X-ray 
luminosity of the cloud, T c is the mean cloud temperature, V c is the mean cloud volume, 
and the subscript represents the initial (preshocked) value. If we simply assume an ideal 
gas so T c oc Pc/pc, we have L x /L xfi ~ (p c / Pc,o) 3/2 {Pc/p c ,o) 1/2 (V c /V Cj0 ) , allowing us to estimate 
the total X-ray luminosity of the cloud with respect to its preshocked X-ray luminosity. In 
model M4, for example, p c / p c ,o ~ 3, p c /p c ,o ~ 12, and V c /V Cj0 « 1 at t/t hc = 3 (as can 
be roughly estimated from Figs. 4 and 5), so that L x /L Xj0 ~ 18. Thus, the shocked cloud 
could possibly be a important source of thermal X-rays depending upon the various physical 
parameters such as the incident jet velocity, the cloud density, and, most importantly, the 
initial cloud temperature, which is not explicitly specified in our scaled models. However 
any thermal significant X-ray luminosity should subside rapidly after the interaction as the 
cloud is diffusive and quickly attains equilibrium with the postshocked ambient pressure. 



5. Summary and Discussion 

We have performed three-dimensional relativistic hydrodynamic simulations to study 
relativistic jet interactions with dense clouds, focusing on the influence of special relativistic 
effects. We have investigated clouds struck by both low and high relativistic beam Mach 
number jets which have less and more dominant relativistic effects, respectively, and have 
compared our results to the extent possible with nonrelativistic simulations which have been 
published previously. We also have studied the evolution of the assaulted clouds and have 
estimated the synchrotron emission from the relativistic jets interacting with the clouds. 

In our models, the partial deflections of the jets due to the interactions with clouds 
are seen more clearly when denser clouds are involved, and the deflection angle is more 
strongly influenced by the density contrast of the cloud to the ambient medium than by the 
beam Mach number of the jet. The streams of deflected gas from the jet induced by the 
interactions move outward much faster compared to nonrelativistic models. If our models 
can be generalized, this suggests that the relativistic jet-cloud interactions are an effective 
mechanism of producing at least some of the outflows observed in the vicinity of AGNs (e.g., 
Emonts et al. 2005; Morganti et al. 2005). After the relativistic jets interact with the dense 
clouds, we find that the slower relativistic jets can be bent by modest angles and that these 
bent jets still remain stable and collimated over fairly extended timescales. This trend is 
similar to the results from nonrelativistic simulations. 

The impact of the jet erodes the cloud, but much of the cloud mass survives as a 
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large coherent body rather than being completely destroyed. This enhancement of the cloud 
durability compared to interactions with planar shocks appears to be primarily due to the 
geometric influence of an off-axis collision. Compared to head-on collisions, off-axis collisions 
damage the cloud less, increasing the chance of survival of a large portion of the cloud. 
Another likely reason for the enhancement of the cloud durability is that the rate of cloud 
fragmentation through Kelvin-Helmholtz instabilities is lowered since the relativistic flows 
reduce the growth rate of the instabilities compared to similar off-axis blows by nonrelativistic 
jets. This leftover tenacious cloud material could be a candidate for strong star formation 
region in the vicinity of AGNs, particularly when the cooling timescales are sufficiently short. 

The synchrotron intensity "maps" show that at the jet impact on a cloud, the syn- 
chrotron emission comes dominantly from a bright hotspot which could correspond to the 
form of the compact emission knots seen in many VLBI radio maps. Although the emission 
from the deflected jet material is relatively weak, there is a secondary synchrotron spot vis- 
ible from this deflected material. This emission feature may represent some of the distorted 
emission seen in many VLBI radio maps. The passage of the jet over a cloud significantly 
enhances the total synchrotron intensity of the jet. We find that the synchrotron emission 
is steeply enhanced shortly after the jet hits the cloud, but the emission peaks right before 
the jet passes through the cloud. The next big step in performing these calculations would 
be to include magnetic fields. Such relativistic magnetohydrodynamical simulations would 
allow for better estimates of synchrotron emission and would allow the examination of the 
polarization of emission arising from more complicated shock structures. Such polarization 
structures could be a useful diagnostic of the dynamics. 

Although most astrophysical simulations based on relativistic hydrodynamics, including 
this study, have assumed the ideal EOS, it is well known that the ideal EOS is correct only if 
the gas is assumed to be entirely nonrelativistic (7 = 5/3) or ultrarelativistic (7 = 4/3). If a 
local transition between nonrelativistic gas and relativistic gas is involved, the ideal EOS will 
produce incorrect results in that regime. Recently, Ryu et al. (2006) have studied this issue 
of the EOS in numerical relativistic hydrodynamics and propose a new EOS which is simple 
and yet approximates closely the EOS of a perfect gas in the relativistic regime, having an 
accuracy in enthalpy better than 0.8%. Future numerical simulations using this new EOS 
should produce even better results concerning the problem of relativistic jet interactions with 
clouds. 

We thank the anonymous referee for several suggestions which improved the presentation 
in this paper. The simulations were performed on the Linux Biocluster at Georgia State 
University, and we are grateful for the allocation of substantial time on this cluster. EC and 
PJW have been supported in part by a subcontract to GSU from NSF grant AST 0507529 
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Table 1. Simulation parameters 



Model 


X 


V 


Vb 


r 6 






t hc 


^end 


Ml 


10 


0.1 


0.9 


2.29 


2.92 


6.36 


4.86 


4t bc 


M2 


100 


0.1 


0.9 


2.29 


2.92 


6.36 


4.86 


6tbc 


M3 


10 


0.1 


0.99 


7.09 


1.92 


11.6 


2.50 


4t bc 


M4 


100 


0.1 


0.99 


7.09 


1.92 


11.6 


2.50 


5tbc 



Note. - Here x is the ratio of the cloud density to the 
ambient medium density, rj is the ratio of the beam density 
to the ambient medium density, Vb is the initial beam velocity, 
r& is the beam Lorentz factor, is the Newtonian beam 
Mach number, M.^ is the relativistic beam Mach number, % c 
is the beam crossing time, and t cn d is the time at which the 
simulation is ended. 
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Fig. 1. — (a) Gray-scale images of density, pressure, and Lorentz factor (top to bottom) for 
model Ml at three different times, t/% c = 1, 2.5, and 4 (left to right). The image scales are 
logarithmic for p and p but the square-root of the logarithm for T so as to enhance visibility 
of intermediate values; the images show the x — y plane with z = 4 in the three-dimensional 
computational domain, (b) Same as Fig. 1(a) except for model M2 at t/% c = 1, 3.5, and 6. 
(c) Same as Fig. 1(a) except for model M3 at t/tb c = 1, 2.5, and 4. (d) Same as Fig. 1(a) 
except for model M4 at t/tbc = 1, 3, and 5. 
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Fig. 2. — Distributions of density, pressure, and Lorentz factor along the beam propagation 
axis with y = z = 4 for models Ml (left) and M3 (right) at the same epoch as in Figs. 1(a) 
and (c). The solid lines correspond to t/t\, c = 1, dotted lines represent profiles at t/% c = 2.5, 
and dashed lines illustrate quantities at t/tbc = 4. 
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Fig. 3. — Projection images of Lorentz factors at viewing angles of 0° for models M2 (top) 
and M4 (bottom) at t/t hc = 3.5 (M2) and 3 (M4). The image scales are logarithmic and the 
images are projected on the y — z plane in the three-dimensional computational domain. 
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Fig. 4. — Volume-rendering images of cloud density for model M4 at three different times, 
t/thc — 1; 3, and 5. The image scales are linear and the viewing area is rotated 20° about 
the x-axis and 30° about the z-axis. Black represents the lowest values which are ~ at 
each epoch and white the highest values which are ~ 100 at t/t\, c — 1, ~ 464 at t/t\, c = 3, 
and ~ 229 at t/t hc = 5. 
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Fig. 5. — Time evolutions of the root mean square radius of the cloud, the mean cloud 
velocity, and the mean thermal energy of the cloud for models Ml (curves ending at t/% c = 4) 
and M2 (curves ending at t/% c = 6) in the left panels and for models M3 (curves ending at 
t/tbc — 4) and M4 (curves ending at t/% c = 5) in the right panels. The solid, dotted, and 
dashed lines in the root mean square radius and the mean velocity panels represent the x-, 
y-, and ^-components, respectively. 



Fig. 6. — Contour maps of synchrotron intensity at the viewing angles of 90°, 45°, and 0° 
(top to bottom) for models Ml (left) and M2 (right) at t/t hc = 2.5 (Ml) and 3.5 (M2). 
Maximum synchrotron intensities are 0.17 (90°), 0.13 (45°), and 0.27 (0°) for model Ml and 
0.57 (90°), 0.35 (45°), and 0.55 (0°) for model M2, and the contour levels are 0.1%, 0.5%, 
1%, 3%, 6%, 10%, 20%, 40%, 70%, and 90% of the maximum synchrotron intensity. 




Fig. 7. — Total synchrotron intensity curves for models Ml (top left), M2 (top right), M3 
(bottom left), and M4 (bottom right). The solid, dotted, and dashed lines correspond to the 
viewing angles of 90°, 45°, and 0°, respectively. 



